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Abstract 

The dual approach to the Ginzburg-Landau theory of a Bardeen-Cooper- 
Schrieffer superconductor is reviewed. The dual theory describes a grand 
canonical ensemble of fluctuating closed magnetic vortices, of arbitrary length 
and shape, which interact with a massive vector field representing the local 
magnetic induction. When the critical temperature is approached from below, 
the magnetic vortices proliferate. This is signaled by the disorder field, which 
describes the loop gas, developing a non-zero expectation value in the normal 
conducting phase. It thereby breaks a global U(l) symmetry. The ensuing 
Goldstone field is the magnetic scalar potential. The superconducting-to- 
normal phase transition is studied by applying renormalization group theory 
to the dual formulation. In the regime of a second-order transition, the critical 
exponents are given by those of a superfluid with a reversed temperature axis. 



I. INTRODUCTION 

In a recent letter [JTJ] critical indices of the superconducting-to-normal phase transition 
were derived starting from a dual formulation of the Ginzburg-Landau theory. The purpose 
of the present article is to review this approach. The article, which is pedagogical in tone 
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and attempts to be self-contained, highlights also some recent developments related to the 
symmetry aspects of the dual formulation which are important for understanding the phase 
transition. 

The basic idea of the dual approach originates from three-dimensional (3-D) lattice studies 
carried out more than 15 years ago These studies were instigated by the success of the 

Berezinskii-Kosterlitz-Thouless (BKT) theory of the 2D xy model developed in the early 
seventies to describe the phase transition in a superfluid film ||[7]]. It was argued that the 
phase transition has its origin in the behavior of point-like vortices. In the ordered low- 
temperature phase, pairs of vortices with opposite topological charge are tightly bound. 
Above the critical temperature, these pairs unbind and the liberated vortices disorder the 
system, turning the superfluid into a normal fluid. 

The lattice studies of the pure 3D xy model and its extension to a lattice supercon- 
ductor |||4j] were aimed at obtaining an analogous description of the phase transitions in 
these models in terms of vortex loops. A detailed account of these matters as well as an 
extensive list of references to the literature can be found in the textbook by Kleinert ||. 

Another development illuminating the dual approach to the Ginzburg-Landau model 
was initiated in Ref. The basic observation was that since a local gauge symmetry 
can never be broken [|lOj, a local gauge description of a phase transition is not feasible. It 



was argued that the 3D Ginzburg-Landau theory contains, in addition to the local gauge 
symmetry, another global U(l) symmetry. When considered in 2 + 1 space-time dimensions 
with a Minkowski metric, this symmetry is generated by the magnetic flux operator. It 
was demonstrated that this symmetry is broken in the normal phase, while it is unbroken 
in the superconducting phase. An order parameter was given, and it was shown that the 
massless photon of the normal phase is the Goldstone particle associated with the broken 
flux symmetry. 

To clarify the concept of duality it is expedient to remain for the moment in 2 + 1 space- 
time dimensions. There the theory is a quantum field theory possessing particle states 
which are created by field operators. This is to be contrasted with the space-time approach 



of Feynman JTTJ] where worldlines of particles are the fundamental objects. The particle's 
mass is given by the tension of the worldline. Now, suppose that a quantum field theory 
contains point-like solitons, i.e., topologically stable, time-independent solutions of the field 
equations. When set into motion such a soliton traces out a non-trivial worldline in space- 
time. An operator description of these excitations in terms of the original fields is non-local 
and quite difficult. The aim of the dual approach of the type we are considering is to describe 
them by a local field, called a disorder field |§. The space-time description of the particle 
is thus converted into an equivalent field-theoretic one. 

In the context of the (three dimensional Euclidean) Ginzburg-Landau theory the notion 
of a particle state is of course missing. The above picture has, however, a direct counterpart, 
viz., worldlines of the space-time approach correspond to vortices in the Euclidean theory. 
More precisely, in the setting of statistical (as opposed to quantum) field theory, a dual 
theory describes vortices by a disorder field. (In this discussion we restricted the notion of 
duality to line defects, thereby omitting dual theories of other topological defects such as 
domain walls and — in the context of quantum field theory — instantons.) 

The advantage of the dual theory is that it allows for a simple description of the su- 
perconducting phase transition in the sense of Landau which is impossible in the original 
formulation, as we shall discuss. The order parameter of the transition is the field repre- 
senting the magnetic vortices. This Landau theory, involving only a global symmetry and 
a genuine order parameter, can — as usual — be taken as a convenient starting point of a 
renormalization group analysis of the superconducting phase transition JT|. 

The layout of the paper is as follows. In the next section we illustrate the dual approach 
by considering the Ginzburg-Landau model in two space dimensions. This example gives 
us important clues towards the form of the dual formulation in 3D. In Sec. [TTT] we review 
the fact that a loop gas can be succinctly described near a critical point by a \tp\ 4 field 
theory. This equivalence plays a central role in the subsequent part of the paper. In Sec. 
P"V| we discuss the dual formulation of a 3D lattice superconductor, and in Sec. [V] we review 



the formulation in the continuum. In Sec. [VJ we argue that the dual theory yields the 



standard description of the mixed state of a type-II superconductor. As a consistency check 
we consider in Sec. |V11| the dual formulation of the dual theory and demonstrate that the 
resulting theory is the original Ginzburg-Landau model. In Sec. |V111| we discuss the gauge- 
invariant superconducting order parameter. Finally, in Sec. [IX] we apply renormalization 
group theory to the dual theory to show that the superconducting phase transition is of 
second order. The critical exponents are shown to be in the universality class of the xy 
model with the temperature axis reversed. 



II. DUAL MAP OF A 2D SUPERCONDUCTOR 

Before discussing the theory of a superconductor in three space dimensions, it is instruc- 
tive to first study the two-dimensional case. This exercise will exhibit the relevant variables 
in a dual description and show what in principle can be expected from such a description. 

The Ginzburg-Landau model resulting from the Bardeen-Cooper-Schrieffer (BCS) theory 
of superconductivity is specified by the Hamiltonian 
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H = d 2 x 



^(V x A) 2 + |(V - 2*eA)0| 2 + mj|0| 2 + A|0| 4 



(1) 



with A the electromagnetic gauge potential. We adopt the standard notation and denote by 
V x A = tijdiAj the local induction. The symbol B is reserved for the spatial average of the 
local field, and will be referred to as the magnetic induction, or flux density. These quantities 
have only one component in two space dimensions, pointing in the missing third direction. 
The coupling 2e accounts for the double charge of Cooper pairs. The complex scalar field <fi 
is the superconducting order field with a mass parameter and a self-coupling A. In the 
superconducting phase of the model, < 0. 

A useful limit of the theory is the so-called London limit where m 2 ^ —>■ — oo and <fi can 
be written as <p(x) = exp[i9(x)]w /a/2, with a constant w. This limit properly represents 
many relevant properties of the superconducting state. The functional integral over the size 
fluctuations of the scalar field can in this limit be approximated by the saddle point, and H 
becomes, omitting irrelevant terms 



T P s 2 + -(Vx A) 2 



(2) 



with p s the superfluid momentum 

p s = V0 - 2eA + P . (3) 

A so-called plastic field 6 P must be introduced for the following reason. The variable 9 is 
a phase variable and is thus cyclic with periodicity 2tt. A jump in 9 by 2ir along a cut 
in the spatial plane indicates the presence of a point vortex. When this is the case the 
Hamiltonian requires a field 6 P to compensate for these jumps so that H remains smooth. 
Such a construction with a compensating field is well known in the Villain formulation of 
lattice models, and it can also be given in the continuum ||. The curl of 6 P gives the density 
of vortices consisting of 5 functions at each vortex position 

V x 6 P = 27r^n Q 5(x-x a ), (4) 

a 

where n a is the winding number of a vortex, and x a its location. We will restrict ourselves 
to vortices with unit winding number, n a = ±1 for a vortex and antivortex, respectively. 
Since the theory is charged, vortices carry also a magnetic flux quantum n a 7r/e, and are 
more properly called magnetic vortices. 
The solution of (U) is given by |§ 

Of = -2vr n ^ij <5j( x > L a ). (5) 

a 

where 5j(x, L a ) is the 5 function on the line L a which starts at the position of the ath vortex 
x a and runs to spatial infinity along an arbitrary path: 

5i(x,L a )= dyiS(x-y). (6) 

With the help of Stokes' theorem in two dimensions which can be expressed in differential 
form as V • <5(x, L a ) = <5(x — x a ) it is easy to check that the explicit expression (^) for the 
field 9f indeed solves Eq. (|J). The shape of the line L a is physically irrelevant, only the 



starting point which marks the position of the vortex matters. Leaving this point fixed the 
shape of L a may be changed at will by performing the transformation 

C(x)^f(x) + «^(x,S), (7) 
where <5(x, S) is the 5 function on the surface S swept out by L a : 

5(x, S) := / dadr ey d a yid T yj S[x - y(a, r)], (8) 
J s 

with y(cr, r) being a parameterization of S. The superfluid momentum and thus the Hamil- 
tonian is invariant under these so-called vortex gauge transformations ||. 

For the time being vortices are ignored by setting the plastic field 6 P to zero. The 
partition function of the system is then given by 

Z = fve J VAE(A) e~ H , (9) 

with H(A) a gauge-fixing factor for the gauge field A, and H the Hamiltonian (0). Here, 
fields and coupling constants are rescaled so that no explicit temperature dependence appears 
in the Boltzmann factor. In ([5]) it is easy to integrate out the 9 fluctuations. This yields for 
the partition function 



Z = J PAS(A)exp|-i| d 2 x 



(V x A) 2 + m\M L - ^) A 3 



(10) 



where the last term with itla = 2ew is a gauge- invariant, albeit non-local mass term for 
the gauge field generated by the Higgs mechanism. A massless gauge field in 2D represents 
no physical degrees of freedom. In Minkowski space-time, this is easily understood by 
recognizing that in 1 + 1 dimensions there is no transverse direction. The system contains 
therefore only a single physical degree of freedom before the Higgs mechanism took place, 
namely 6. This equals the number afterwards since a massive vector field represents only 
one independent degree of freedom in 2D. Note that the absence of genuine long-range order 
in 2D is no obstacle to the Higgs mechanism. (The massless Schwinger model ]T2|] provides 
the simplest soluble example for this.) 
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We next introduce an auxiliary field h to linearize the first term in fllQ|) 

ir„,- . ,oi /■_. r i 



cxp 



-- J d 2 x(V x A) 2 = J Vh exp -- J d 2 x h 2 + i J d 2 x h • ( V x A) 



and integrate out the gauge-field fluctuations [after adding a gauge-fixing term (1/2£)(V • 
A) 2 ]. The result is a manifestly gauge-invariant expression for the partition function in 
terms of a massive scalar field h, which represents the single degree of freedom contained in 
the theory: 



£>hexp J -- / <P.r 



^(Vh) 2 + h 2 



(12) 



To understand the physical significance of the h field appearing in this functional integral, 
we note that it follows from ([TTj) that the auxiliary field h satisfies the equation 



iV x A. 



(13) 



That is, the fluctuating field h coincides up to a factor i with the local induction. Equation 
flT2|) shows that the magnetic field has a finite penetration depth A = \jm^. In contrast to 
the Ginzburg-Landau description where the functional integral runs over the gauge potential, 
the integration variable in ([12]) is the physical field. 

We now include vortex contributions. The mass tjia provides the system with an infrared 
cutoff so that a single vortex in the charged theory has a finite energy, implying that there 
will always exist thermally activated vortices. This is different from the neutral model, 
describing a 2D superfluid, where the absence of an infrared cutoff permits only tightly 
bound vortex- antivortex pairs to exist. We expect, accordingly, the superconducting phase 
to describe a plasma of vortices, each carrying one magnetic flux quantum ±n/e. The 
partition function now reads 
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E 



J] d 2 x a I V6 I VA exp \ ~ I d' z x [w 2 p z s + (V x A) 



(14) 



n + 7Z=0 n + lU - 

where the superfluid momentum p s contains the vortex gauge field P , as in (^]). In 
(HD> IIlad 2 x a denotes the integration over the positions of the vortices and the factor 



l/n+\ (l/n_!) accounts for the fact that the (anti-)vortices are indistinguishable. The fu- 
gacity z is the Boltzmann factor for an isolated vortex z = exp(— e c ), with e c being the 
vortex core energy. 

The field 6 P can be shifted from the first to the second term in ( |T3D by applying the 
transformation A — > A — p /(2e). This results in the shift 



VxA^VxA-5 1 



(15) 



with the plastic field 



x — x r 



(16) 



representing the magnetic flux density. Repeating the steps of the previous paragraph we 
now obtain instead of (|12|) 



E 



n + ,n_=o n+\n_\ 



Phexp { -- / d 2 x 



^(Vh) 2 + lV 



+ i / d 2 xB p h\. (17) 



where h represents the physical local field 



h = z(V xA-B 1 



(18) 



as will be clarified in a later section (Sec. [Vj). The last term in fljj]) shows that the charge 
g with which a magnetic vortex couples to the fluctuating h field is the product of an 
elementary flux quantum n/e (which is contained in the definition of B p ) and the inverse 
penetration depth mi, i.e., 



7T 

g = -m A . 

e 



(19) 



For small fugacities the summation indices n + and n_ can be restricted to the values 0, 1 



and we arrive at the partition function of the massive sine-Gordon model [13 



Vh cxp 



d x 



1 



(Vh) 2 + h 2 



"2.X cos ( -h 



(20) 



This is the dual formulation of a 2D superconductor. The vortices of unit winding number 
n a — ±1 turned the otherwise free theory fll2|) into an interacting one. 
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The final form ( p0|) demonstrates the basic concepts of a dual formulation: it is a formu- 
lation directly in terms of a (gauge- invariant) field representing the physical local induction. 
This is different from the Ginzburg-Landau description of a two-dimensional superconductor 
where the magnetic field is the curl of an unphysical gauge potential A. The dual formu- 
lation also accounts for the topological excitations, viz., the magnetic vortices which in 2D 
are point defects. They are coupled with a charge g = (ir/e)mA to the magnetic field. 

In three space dimensions these two basic ingredients will surface again. Since in 3D the 
magnetic field has three components of which two are independent, the dual formulation will 
involve a massive vector field, rather than a simple massive scalar field h as in 2D. But what 
is more important, the point vortices of the two dimensional case become line defects in 3D. 
A grand canonical ensemble of closed loops requires, as will be demonstrated in the next 
section, a fluctuating field of its own to specify it. That is, accounting for the vortices does 
not simply add an interaction term to the theory without vortices, but adds a whole new 
theory. The coupling of this new theory — which turns out to be a complex |?/>| 4 theory — to 
the field representing the local field is again given by (7r/e)mA as in two dimensions. 



III. FUNCTIONAL-INTEGRAL DESCRIPTION OF A LOOP GAS 

The subject of this section is the well-known fact that a loop gas can be described with 
the help of a functional integral involving a complex theory P,|H|. This will be used in a 
later section when formulating the dual theory of the Ginzburg-Landau theory in 3D. This 
theory features fluctuating closed magnetic vortices described by a disorder field theory. We 
will review here the derivation in the continuum, the discussion of the lattice derivation is 



relegated to the Appendix |8|,|15H 



Our point of departure is the correlation function G(x) of a free complex field theory 
defined by the Hamiltonian 

H = J d 3 x(|VVf + mV| 2 ) • (21) 

Explicitly, 



d 3 k c 



,'ik-x 

(2vr) 3 k 2 T^i 2 



^(x) = / ( 22 ) 



This is written in the Schwinger proper-time representation as an integral over the proper 
time s |nj 



' '' d3k ,k,. s k* 



2 / 1 1 2/ 



cfce— ^— ) e~r", (23) 
where we used the identity 

- = / dse~ sa . (24) 
a Jo 

According to Feynman's formulation of quantum mechanics |TJ] the right-hand side of ( |23|) 
can be represented as a sum over all real-space paths of a particle with mass \ running from 
at imaginary time to x at time s: 

(jM e-I x2 /^= / Dx( S ')exp _i f rfs 'x 2 ( S ') , (25) 

V47TS/ Jx(0)=0 L JO 

involving the free "Lagrangian" L = x 2 /4. The extra Boltzmann factor exp(— sm 2 ) in ( |2"3"| ) 
suppresses large proper-time values exponentially, and the integral J °° ds shows that all 
positive values of s are allowed and accounted for. 

From integrating the identity (|4]) with respect to a it follows that to within an additive 
constant 

ln(a) = - / — e~ sa . (26) 
Jo s 

Employing this we can use the previous result to write 



Tr ln(-V 2 + m 2 ) = - f°° — e"™ 2 / -^-e 



ds_ Q - sm i f d 3 k 
os J (2vr) 

00 ds 
o s 



= 'Jo" s " m2 / Px(s,)eXP -lJ Sds ' ±2 ^ ' (27) 
where Tr denotes the integral over the momentum variables 
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Trln(-V 2 + m 2 ) = / -f^- ln(k 2 + m 2 ). 

J yllt) 6 



(28) 



The path integral § Ux(s') in ( P7|) runs over closed loops starting and ending in 0. Using 
the identity 



Trln(-V 2 + m 2 ) = In Det(-V 2 + m 2 ), 



(29) 



we can write the inverse determinant as 



Det 



-V 2 + m 2 ) = exp(-Wo) 



(30) 



with Wo denoting the right-hand side of (|27j). The exponential exp(— Wo) has the expansion 



-Wo 



oo I N 
N=0 iv • 1=1 



exp 



N 



^ z=i ,yu 



(31) 



/0 S; 7 

which we recognize as the partition function Z of a grand canonical ensemble of fluctuating 
closed oriented random loops, of arbitrary length and shape. We will refer to such an 
ensemble as loop gas. On account of (^), the exponential exp(— Wo) may alternatively be 
viewed as an inverse functional determinant. As such it is easily recognized as the partition 
function of the free complex field theory, 



Z = J V^*Vi)e~ H , 



(32) 



with H the Hamiltonian Thus we have proved the equivalence of a free complex field 
theory and a free loop gas. 

In a superconductor, the vortices have short-range interactions. To describe these we 
start again from the field theoretic side. We shall prove that the Hamiltonian 



H 



d 3 x 



|W| 2 + m>| 2 + A(H 



1 2\ 2 



(33) 



with the additional interaction term, accounts for the steric repulsion in the loop gas. For 
this we write the complex field ip as tp = (ipi + iip2) / V% and express the interaction term as 
a functional integral over an auxiliary field a 
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exp 



2\2 



"Dcxexp 



d 3 x Qa 2 - # 2 o) 



(34) 



where the index a runs over 1,2. The partition function becomes now 
Z = I Vc.Vc-.Vavx]) { - I d 3 x 



J DcrDet _1 (-V 2 + m 2 - 2ia) exp (~jf d 3 xa 

I v, exp (-i / ,>) t m n [f f 



X 



{ N r*i rl 



(35) 



where the last equality follows from the previous result ([H]). A simple Gaussian integration 
yields 

J Vaexpl^-jJ (fW + 2!^s'ff[x(s')]} = exp |-A J" ds'ds" 8 [x(s') - x(s")]} . (36) 

Using this in the last expression for Z, we obtain the real-space representation for the 
partition function of a complex |?/>| 4 theory 



N 



oo 1 
jV=0 iv • i=l 



iV 



exp ( _ 4 ^ Jo rfs ' i2 ^'' 



-A E f ^ r ^ [x(*0 - x(4')]| • (37) 

l,k=l J ° J ° J 



This is recognized as the partition function of a loop gas with short-range repulsion. The 
field theoretic representation of an interacting loop gas by a |?/>| 4 field theory will be employed 
in Sec. |V| to formulate the dual theory of a 3D superconductor. 



IV. DUAL TRANSFORMATION OF A 3D LATTICE SUPERCONDUCTOR 

In this section we set up a dual formulation of a 3D superconductor on a lattice. In 
contrast to the continuum case to be discussed in the next section, the lattice model can be 
transformed exactly to its dual version. We take as starting point the partition function |3j 
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rfA(x)5(V ■ A) 



e -IEx(vxA)2-Q r r 



lw/a 2n 



exp < (3 cos[<9j#(x) 



gA(x)] 

(38) 



where 8(x) is a phase variable at site x of a 3D cubic lattice, A(x) is the electromagnetic 
gauge potential which is a real non-compact variable defined on the directed links between 
adjacent sites, (3 is the inverse temperature, and q is the electric charge. Moreover, x = axii, 
with a the lattice spacing, x% integers labeling the sites, and i three orthogonal vectors 
spanning the lattice. The sum J2 X , which includes a factor a 3 , runs over all lattice sites, 
while the sum J2i runs over all directions; V is a lattice derivative with components 



$/(x) = ~[/(x + oi) - /(x)]; 0,/(x) = i[/(x) - /(x - ai)], 



(39) 



and the delta function 8( V • A) in (p8|) fixes the gauge. The lattice model (|38|) is appropriate 



for a non-compact gauge group. We shall be working in the Villain approximation |17] of 
the model, which is obtained by the following replacement 



n 

X 

n 



rn/a 


dd(*)~ 


Jrc/a 




rn fa 




J tv /a 


2tt 



exp 1(3 cos[^6»(x) - gAj(x)] 

^ i,x 

E exp{-§5:[V^(x)-gA(x)-27rn( 



*)f 



(40) 



where n(x) are integers which are — like the gauge potential A(x) — defined on the directed 
links between adjacent sites. The partition function in (|40"D is given the index xy to indicate 
that this part of the theory is simply the xy model when the electric charge q is set to zero. 
We proceed by introducing an auxiliary field v(x) via a quadratic completion 

^ d0(x) 



n 



ir/a 2n 



dv(x) 



x 



E ex P {-^E v 2 (x) + ^E v(x) ■ [W(x) - gA(x) - 2vrn(x)]] . (41) 

{n(x)} ^ ' x x / 



The integration over 9{x) can then be carried out to yield the constraint V • v(x) = 0, while 
the sum X){ n (x)} exp[— 27riv(x) • n(x)] forces the integral over the real variable v to take on 
only integer values 1, say. This last observation follows from the Poisson summation formula 
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-2nivn 



(42) 



In this way the partition function Z xy becomes: 



Z xy = E %l,O eX P 
{l(x)} 



-^E l2 W-^E!(x)'A(x) 

Z P X X 



(43) 



Whereas representation fliOD of Z xy is a field theoretic description of the lattice model, Eq. 



(f43|) describes it in the real-space language of (charged) closed loops. This is analogous to 
the twofold description of a loop gas in the continuum which we discussed in the preceding 
section. 

The constraint V • 1 = in fl4*3[ ) can be explicitly solved by introducing an auxiliary 
integer- valued potential i(x), via 1 = V x i, so that 



n 



dh(x)5(V- h) 



E % m , e- 2 ^ m - h exp 

{m(x)} 



— X:(Vxh) 2 -zgX:(Vxh).A 

(44) 



Here, because of the presence of the factor exp(— 2ni Y^ x m • h), the sum over the integer- 
valued fields m(x) ensures, on account of the Poisson formula fl42|) , that only integer- valued 
fields h(x) = i(x) contribute to the integral / dh(x), as required. For q — 0, corresponding to 
the uncharged xy model, Eq. (f44|) gives a dual representation in terms of vortices m coupled 



to a massless vector field h. The condition V • m = in (|il) is to assure that the vortices 
form closed loops, which is physically the case. When the integration over the field h is 
carried out, so that one is only left with vortices, one finds them interacting via long-range 
forces of the Biot-Savart type. 

For the lattice superconductor we have to add the electromagnetic part to Z xy [|J: 



n 

X 

x exp 

n 



dA(x)8(V ■ A) / dh(x)5(V- h) 



E ^V-m.O 
{m(x)} 

-\ E(V x A) 2 - 27TZ E m • h - i E(V x h) 2 - %q E(V x h) • A 

^ X X X X 



dh(x)5(V- h) 



E %. m , e- 2m ^ m ' h exp^-^E [(V x h) 2 + m 

{m(x)} y ' 



(45) 
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where in the last step we carried out the integral over the electromagnetic gauge potential 
which resulted in a mass term for the vector field h. The constraint Vh = is an intrinsic 
part of the description of a fluctuating massive vector field. For a non-fluctuating field 
it follows automatically from the field equation. Physically, ( |45| ) represents a loop gas of 
vortices coupled to a massive vector field. The massiveness of h is in fact the only difference 
with the dual description of the xy model, Eq. ([y]) with q = 0. As a result, when integrating 
out this vector field, one obtains again a Biot-Savart type of force between the vortices, but 
now of finite range. 

In the next section, an equivalent dual map will be carried out in the continuum, yielding 
a continuum version of the lattice model ([45]) . The fluctuating vector field h will be related to 
the physical local induction. With this identification in mind the constraint V • h = in ([44]) 
can be understood as representing the fact that the magnetic induction is divergence- free. 



V. DUAL MAP OF A 3D SUPERCONDUCTOR 

This section is the central part of the paper, in which we derive the dual formulation 
of the Ginzburg-Landau model. We shall argue that this description is one in terms of a 
genuine order parameter which involves a global rather than a local symmetry as is the case 
in the Ginzburg-Landau formulation ||. For this reason the dual theory can be employed to 
arrive at a conventional Landau description of the superconducting phase transition p|,|l8|,|l9| , 
which will be the subject of Sec. |X[ 

The fundamental object of the dual description is the magnetic vortex, or Abrikosov 
flux tube. Such a topological defect is either closed, infinitely long, or — in the case of a 
finite system — starts and ends at the boundary of the specimen. A magnetic vortex can 
never terminate inside a superconductor. However, for reasons that shall become clear when 
we proceed, we will employ a construct that allows us to describe, at least theoretically, 
magnetic vortices that do terminate inside the superconductor. To this end we allow the 
system to contain a Dirac monopole []2D| as a test particle. Recall that due to the Meissner 
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effect, a magnetic field can penetrate a superconductor only by forming quantized flux tubes. 
Since also the flux lines emanating from the monopole are squeezed into a tube, a monopole 
produces precisely such a vortex. 

Section |I| revealed that a magnetic vortex is described by a plastic field B p which appears 
in the theory in the combination V x A — B p with the gauge field |2t| . In other words, in 
the presence of a test tube, the Hamiltonian becomes 



d 3 x 



i(V x A - B p ) 2 + \rn\M L - 9 -0\ Aj + 1(V • A) 2 



(46) 



where we added a gauge-fixing term (1/2£)(V- A) 2 , and given H the superscript P to indicate 
the presence of the test particle with its emanating tube. We recall that we consider the 
Ginzburg-Landau model (|l|) in the London limit, where the superconducting order field is 
written as <f)(x) = exp[i9(x)]w / 'a/2, with w a constant. The mass term, with = 2ew, is a 
result of integrating out the phase field 9. In three dimensions, the plastic field B p has the 



form 20.21 



fi p (x) = -/ d A[x - y( a )] = - / dy, «5(x-y), 
e jl z as e Jl 2 



(47) 



which is the proper three-dimensional generalization of the 2D result ([H]). The field satisfies 
the equation V • B p (x) = (7r/e)5(x — z) on account of Stokes' theorem, implying that it 
indeed describes a monopole located at z. The line L z starting at the point z and running 
to infinity is the Dirac string accompanying the monopole. As will be clarified below, the 
location of the flux tube coincides with that of the Dirac string (see Fig. |l|). 

From ( f46"D we infer that the operator V(L Z ) describing the test tube is given by 



V(L 2 



e 2 



exp 



J d 3 x (V x A) ■ B 1 



(48) 



We are interested in the expectation value of this operator: 



(V(L.)) = J VAe 



(49) 



Since the integral over A is Gaussian, it can be evaluated by substituting the field equation, 
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A.(x) = I dV^x - y ) [V x B p (y) 



(50) 



back into the exponent. The gauge-field correlation function A™ appearing here is 



d 3 k ( 5ij (kikj)fk 2 kikj\ ik . x 



(27T) 



k 2 + 



+ c- 



k 4 



(51) 



The gauge-dependent longitudinal part of the correlation function does not contribute to 
(P3) since V- (V x B p ) = 0. The expectation value is therefore independent of gauge choice. 

The local induction corresponding to the classical solution ( j50|) in the presence of the 
test tube is 



V x Afx) -B p (x) 



7T 



VA(x -z)-m 2 A J d 3 yA(* - y)B p (y) 



(52) 



where 



Afx) 



d 3 k 



,tkx 



-m A \x\ 



(53) 



(2%) 3 k 2 + m 2 A ' 4tt|x| 

is the Yukawa potential. The term B p in ( |52| ) describes the Dirac string L z , the first term on 
the right-hand side corresponds to the screened Coulomb force generated by the monopole. 
The last term, which is only present when ^ — i.e., when there is a Meissner effect- 
describes the magnetic flux tube. A closer inspection of ( |52"D reveals that the subtraction of 
the Dirac string B p from the field V x A is necessary in order to obtain the physical local 
induction h [20,^1 1. Indeed, if we calculate from the right-hand side of (|52f) the magnetic 
flux through a plane perpendicular to the Dirac string, we find 



J d 2 x t [^A(x -z)-m 2 A J d 3 yA(x - y)Bf(y) 
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e 



(54) 



that precisely one flux quantum pierces the surface in the negative direction (see Fig. |l]). 
Here, d 2 Xi is an element of the surface orthogonal to the Dirac string, and we used Gauss' 
law to rewrite the first term on the left-hand side as 



J d 2 Xi ftA(x -z) = ^Jd 3 x V 2 A(x - 



z . 



(55) 
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Equation ([54]) confirms Dirac's statement that the magnetic flux emanating from a monopole 
must be supplied by an infinitesimally thin string of magnetic dipoles and that in order to 
obtain the true local field of a genuine point monopole, this string has to be subtracted. 
While this string is indeed completely unphysical in the normal phase, it acquires a physical 
relevance in the superconducting phase |22j where it serves as the core of the Abrikosov flux 
tube. 

Substituting the field equation (|5CiD back into the theory, we obtain for the vacuum 
expectation value (V(L Z )) in the London limit 



(V(L,)> = el/ d M BP ) 2 exp |i J d 3 xd 3 y [V x B p (x)], A y (x - y) [V X B p (y)],} , 



(56) 



or 



(V(L Z )) = exp{~ J d 3 xd 3 y [p(x) A(x - y) p(y) + m 2 A Bf{x) A(x - y) B? (y)] }, (57) 

where p(x) = (7r/e)5(x — z) is the monopole density. In deriving this we have omitted terms 
depending only on w. 

It should be noted that the first factor in (SB) diverges 



(58) 



representing the self-interaction of the Dirac string. This term canceled in (p7|) . The first 
term in ([57]) contains for x = y a diverging monopole self-interaction. This divergence is 
irrelevant and can be eliminated by defining a renormalized operator 

r 2 



V I (L z )=V(L z )exp 



7T" 

2^2 



A(0) 



(59) 



The second term in (|5"7 ) is the most important one for our purposes. It represents a Biot- 
Savart interaction between two line elements dx{ and di/i of the magnetic vortex (see Fig. [|). 
It contains an ultraviolet singularity due to fact that in the London limit, where the mass 
|m^| of the superconducting order field is taken to be infinite, the vortices are considered to 
be ideal lines. For a finite mass the magnetic vortices have a typical width of the order of 



the coherence length £ = l/|m^|. This mass therefore provides a natural ultraviolet cutoff 
to the theory. The last term in fl5T|) then takes the form ||2"2| 



m? A ( 7T x 2 



f dxi [ d Vl A(x - y) = -My \L Z \ , (60) 



with \L Z \ the (infinite) length of the flux tube, and 



TT\ 2 m 2 A , {\ml\\ Tiw 2 , /A 2 " 



M ' = UJ ^ ln (^fJ = — ln (?J (61) 

being the free energy per unit length. The combination A/£ defines the dimensionless 
Ginzburg-Landau parameter k, which in the London limit is much larger than 1. 

For a monopole-antimonopole pair, (|60|) amounts to a confining linear potential between 
the monopole and antimonopole in the superconducting phase. Let V*(L Z ) describe an 
antimonopole located at z, with L z being a line running from infinity to z. Collecting all 
terms, we find for such a pair 

(V r (L z )V;(L z )) = exp(-M v \L ZZ \) exp f— j , (62) 

where L zz is the flux tube connecting the monopole at z with the antimonopole at z, and 
\L ZZ \ is its length. 

We remark that the two Dirac strings may initially run to any point at infinity. Due to 
the string tension, they join on the shortest path L zz between the monopoles. The result ( j62|) 
is central to our line of arguments. It shows that the correlation function (V T (L Z )V* (L z )) 



behaves differently in the two phases pj24| . In the superconducting phase, where the gauge 
field is massive, the "confinement" factor dominates and the correlation function decays 
exponentially for distances larger than 1/My\ 

(V V (L Z )V T *(L,)) |L ^°°o. (63) 

This behavior is typical for an operator in a phase without massless excitations. On the other 
hand, in the high-temperature phase, where the gauge field is massless, the confinement fac- 
tor in the correlation function ( |62"D disappears, while the argument of the second exponential 

19 



turns into a pure Coulomb potential. The correlation function remains, consequently, finite 
for large distances: 

(^(L Z )K*(L,)> |L ^°°1. (64) 

By the cluster property of correlation functions this implies that the operator describing 
the test tube develops a vacuum expectation value. This signals a proliferation of magnetic 
vortices. Indeed, according to (^1|) the free energy My per unit length of a vortex vanishes 
at the transition point, where w — > 0. It should be noted that it is the high-temperature 
phase and not the superconducting phase where V T (L Z ) develops an expectation value. 

Before deriving the full-fledged dual theory let us rederive the correlation function fl6"2"|) 
in a way that reveals some aspects of the nature of the dual theory. To this end we linearize 
the functional integral over the gauge field by introducing an auxiliary field h. In the gauge 
V • A = 0, which corresponds to setting ( = 0, we find 

(V r (L z )V;(L- z )) = J DAPhexp ~Jd 3 xh 2 + i |d 3 ih - (V x A-B p ) - y J d 3 xA 



2 



(65) 

where now V • B p (x) = (7r/e)[5(x — z) — <5(x — z)]. Also the divergence V • (V x A — B p ) 
is non-zero only at the location of the artificially introduced monopoles. In the absence of 
monopoles it follows that only the transverse part of h couples to V x A — B p . We therefore 
restrict the integral over the auxiliary field h to the transverse degrees of freedom. This is 
justified by considering the field equation for h following from (j65|) 



h = i(V x A- B p ) = ih. (66) 

It tells us that apart from a factor i, the fluctuating field h may be thought of as representing 
the local induction h, which we know to be divergence-free in the absence of monopoles. 

The integral over the vector potential is again easily carried out by substituting the field 
equation for A, 

A = -^-V x h, (67) 



ni 



2 
A 
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back into (0), with the result 

{V r (L x )V*(L % )) = [ Vh 5(V • h) exp { -\ [ d 3 x -^(V xhf + h 2 - % [ d 3 x h ■ B 1 



(68) 



We have incorporated a 6 function enforcing explicitly the constraint V • h = 0. In the 
present formulation this constraint is an intrinsic part of the description of a fluctuating 
massive vector field. For a non-fluctuating field it is a consequence of the field equation of 
h: 

- —{didj - V 2 8 ij )h j -hi- iBf = 0. (69) 
m A 

Applying di to this equation, we obtain V • h = provided no monopoles are present and 
the mass is non-zero. 

Expression (|68D shows that the test tube described by the plastic field B p couples to the 
fluctuating massive vector field h, with a coupling constant given by g = {ji/e)m A = 2nw 
as in two dimensions. As T approaches the critical temperature from below, w goes to zero, 
and h decouples from the test tube described by B p . After carrying out the integral over h 
in (|68|) we recover the result (|62|) . 

The fact that the magnetic field has a finite penetration depth in the superconducting 
phase is reflected by the mass term of the h field. 

It is interesting to consider the limit m A — > in detail, where the massive vector field 
decouples from the magnetic vortex. This limit yields the constraint V x h = which can 
be solved by setting h = V7. The correlation function (V T (L Z )V* (L^)) then takes the simple 
form 



(V t (L,)V x *(L t )) = J exp [-Iy> x (V 7 ) 2 + * / d 
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(70) 



where p is the monopole density. In the absence of monopoles, the theory reduces to that of 
a free massless mode 7 that may be thought of as representing the magnetic scalar potential. 
This follows from combining the physical interpretation of the vector field h ( p^| ) with the 
equation h = V 7 . Specifically, 
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V7 = i(VxA-B p ). (71) 

Using the definition of the monopole density, p(x) = (n/e)[5(x — z) — <5(x — z)], we see that 
in terms of the field 7 the correlation function reads 

5 (7r/e)i[7(z)-7(z)] 



(V t (L m )v;(Li) = (e [n/eM ^ z >-^) . (72) 

This demonstrates that the operator V T (L Z ) describing the test tube, which was introduced 
in (f48|) in a real-space formulation involving the singular plastic field B p fl47Q , is now repre- 
sented as an ordinary field. Since we are in the normal conducting phase, where V T develops 
a non-zero expectation value, the presence of the phase 7 indicates that this expectation 
value breaks a global U(l) symmetry, with 7 the ensuing Goldstone field. This will be 
further clarified below. 

Equation (ff2[) reveals in addition that in the normal conducting phase the Dirac string 
looses its physical relevance, the right-hand side depending only on the end points z and z, 
not on the line L zl . This fact is also apparent from our starting formula (|57"D , where the last 
term and therefore any reference to L z disappears in the limit mA_ — * 0. It makes no sense 
to talk about magnetic vortices in this phase because they are condensed and do not exist 
as physical excitations. There is also no non-trivial topology to assure their stability. 

We are now in a position to derive the dual theory of a 3D superconductor. This theory 
features a grand-canonical ensemble of fluctuating closed magnetic vortices, of arbitrary 
shape and length, which have a steric repulsion, i.e., a loop gas of magnetic vortices. We 
know that such an ensemble can be described by a disorder field theory, consisting of a 
complex \ip\ 4 theory. On the other hand, our study of a single magnetic vortex revealed that 
it couples with a coupling constant g to the fluctuating vector field h. These two observations 
uniquely determine the dual theory in the London limit as being given by [25,26,8,18,19 



Z = J VhV^Vij) 5(V ■ h) exp (- jfy) (73) 

with 
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/ 



1 



_ (V x h) 2 + -h 2 + |(V - z-h)^| 2 + m%\^ + u\^ , 



(74) 



x 



2ml 



where the ip field is minimally coupled to the vector field h. Equation (|74]) replaces the 
lattice Hamiltonian ( [45|) near the critical point. It is a description of the superconducting 
state in terms of physical variables: the field h describes the local induction, whereas ip 
accounts for the loop gas of magnetic vortices. There are no other physical objects present 
in a superconductor. The dual theory has no local gauge symmetry because the vector 
field h is massive. In fact, the two observations are connected. The presence of a local 
gauge symmetry in a given theory may be looked upon as reflecting a redundancy in the 
description. Since the dual theory is formulated in terms of physical variables, there is no 
redundancy, and thus no local gauge symmetry. 

Although ( |73|) was derived starting from the London limit, it is also relevant near the 
phase transition. The point is that integrating out the size fluctuations of the scalar field <fr 
would only generate higher-order interaction terms and a possible change of the mass and 
interaction parameter and u. But these modifications do not alter the critical behavior 
of the theory. 

The energy My (^1|) appears in the dual theory as a one-loop on-shell mass correction 
stemming from the graph depicted in Fig. |2|, which we now interpret as a Feynman graph. 
The straight and wiggly lines represent the ip and A field correlation functions, respectively. 

A measure for the interaction strength of a massive vector field in 3D is given by the 
dimensionless parameter equal to the square of the coupling constant multiplied by the range 
of the interaction. For the dual theory this factor is g 2 /rriA ~ ^a/c 2 , which is the inverse 
of the strength of the electromagnetic gauge field A in the superconducting phase. This is 
a common feature of theories which are dual to each other. 

Another notable property of the dual theory is that in the limit e — > it changes into a 
local gauge theory ||, 




d 3 x -(Vxh) 2 +|(V-z^h)^| 2 + m 2 |V| 2 + M |^| 4 , 



(75) 



as can be checked by rescaling the dual field h in the Hamiltonian ([74]) . 
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We next investigate what happens with the dual theory when we approach the critical 
temperature. Remember that w and therefore iua tends to zero in the limit where T tends 
to the critical temperature from below. From the first term in the Hamiltonian ( |74|) it again 
follows that V x h — > in this limit, so that we can write once more h = V7, and flT4"|) 
becomes 



H, h = / d 3 x 



I ( v 7 ) 2 + |(v - i-v 7 )Vf + m 2 M 2 + u|Vf 



(76) 



This equation shows that 7, representing the magnetic scalar potential, cannot be distin- 
guished from the phase of the disorder field. Indeed, let (ir/e)'d be this phase. Then, the 
canonical transformation $ -—>■ $ + 7 absorbs the scalar potential into the phase of ip; the 
first term in ( ff6| ) decouples from the theory and yields a trivial contribution to the partition 
function. In this way, the dual theory reduces to a pure \ip\ 4 theory 

H+ = f £x (|Wf + mjH 2 + u\^) . (77) 

It was already concluded that in the high-temperature phase the magnetic vortices proliferate 
as indicated by the fact that ip, giving a field theoretic description of the loop gas of these 
objects, develops a non-zero expectation value at the transition point. This transition is 
triggered by a change in sign of m 2 . In the London limit the Hamiltonian ( |7T|) then takes 
the simple form 



H^p = I (fx 



(78) 



2 

with v the expectation value of the disorder field, v/\r2 = (\ip\), and where we now repre- 
sented the phase of if) by (ir/e)^ to bring out the fact that 7 describes the magnetic scalar 
potential. As we will demonstrate below, v has the value v = e/ir of an inverse flux 
quantum, so that with our normalization choice of the phase of the if) field, Eq. ( [78]) takes 
the canonical form. 

The picture of the superconducting phase transition that emerges in the dual formulation 
of the Ginzburg-Landau theory is the following. When the critical temperature is approached 
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from below, there is a proliferation of magnetic vortices. We recall that in the London limit 
parallel vortices repel each other, so that a single vortex prefers to crumple. Near T c we 
then have a spaghetti of vortices which fill the space completely at and above the transition 
temperature. Since inside the core of a vortex one has the normal phase, the system thus 
becomes normal conducting. Whereas in the Ginzburg-Landau formulation a magnetic 
vortex is described by a singular plastic field B p , in the dual formulation it is represented by 
the Noether current ji = ijj*diip — 2ighiij)*ip. This follows from comparing the terms coupling 
linearly to the fluctuating h field. In the normal conducting phase the field ip develops a 
vacuum expectation value, and thereby breaks the global U(l) symmetry of the {ipl 4 theory; 
7 is the ensuing Goldstone field. The Noether current becomes in the London limit j = V7, 
with 7 representing the massless photon of the high-temperature phase. It should be noted 
that at T c the fluctuating local field h decouples from ip because g — > 0. 

VI. THE MIXED STATE 

In the previous section it was shown that the dual theory of the Ginzburg-Landau model 
features a loop gas of magnetic vortices, i.e., closed random vortices of arbitrary shape 
and length which are generated by fluctuations. It is well known that magnetic vortices 
can also be generated by applying an external magnetic field H. Above a certain critical 
value H Cl magnetic vortices start to penetrate the superconductor, provided it is a type-II 
superconductor. Below this lower critical field, the Meissner effect expels all flux lines from 
the system. In thermodynamic equilibrium the vortices in the so-called mixed state form 
a 2D triangular lattice perpendicular to the applied field, a so-called Abrikosov flux lattice 
p3fl (see Fig. |3|). The 2D vortex density in such a lattice is given by = 2/ (v^a 2 ), with a 
the lattice spacing. The magnetic flux $ = / d 2 xB through the lattice, which is given by the 
number of vortices N multiplied by the fundamental flux unit 7r/e carried by a single flux 
tube, increases with the applied field. At a second critical value H C2 the magnetic induction 
B becomes homogeneous and saturates the applied field; no more vortices are nucleated. 
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The maximum vortex density p®, m ax reads 



6 i 

P®,max H C2 — — , (^9) 

7T 27T£; 



where H C2 = 1/ (2e£ 2 ) is the upper critical field expressed in terms of the coherence length 
£ = l/lm^l p?fl . At this value the magnetic vortices are closely packed and the system 
becomes normal conducting (see below). The area that can be assigned to a single vortex 
is S® = 1/p® = 27r£ 2 , so that NS® covers the whole surface S perpendicular to the applied 
field. 

Let us study the Abrikosov flux lattice in the London limit, where the thickness of a 
magnetic vortex is considered to be infinitesimal small as compared to the penetration depth 
A = 1/rrijy. In this limit, where |m^| = l/£ may be taken to be infinite, the superconducting 
order field <fi is frozen in, so that there is only one type of interaction between the magnetic 
vortices mediated by the gauge field A. This magnetic interaction is repulsive for two 
parallel vortices. (Outside the London limit, where the coherence length is finite, there is a 
second interaction mediated by <fi. This interaction is attractive for parallel vortices.) The 
Abrikosov flux lattice is easily understood in the London limit. Due to the repulsive magnetic 
interaction the vortices are driven apart, but this is called to an halt by the finiteness of the 
system. The vortices then order themselves in a regular lattice — which turns out to be a 
triangular lattice — so as to minimize the repulsive interaction. 

The lattice is described by the plastic field B cxt which has only a component in say the 
negative third direction: 

5 ext (x) = - - I dx 3 <f(x - Xa), (80) 
a e J 

where the sum is over all lattice points x Q , and where we took the external field as pointing 
in the third direction. In a first approximation, the ip field in ([74]) may be neglected. The 
field equation for the local field h = (0,0, h(x)) which we infer from ( f74|) augmented with 
the term — ih • B ext then reads 

- X 2 V 2 h + h = -B ext , (81) 
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where we accounted for the factor % between the fluctuating field h and the physical local 
field h. The form flSlj ) is well-known, and can also be obtained from the Ginzburg-Landau 
theory. We consider the lattice in two different limits. The first limit is the one of low vortex 
density. Just above the lower critical field H Cl the density is such that one may consider the 
system as non-interacting. In this case the field equation (^T|) has the solution 

Mx) = ^E^o(A- 1 |x-x Q |), (82) 

where K is a modified Bessel function. The corresponding free energy (density) is 

F = -— I d 3 xhB ext 



2V 
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E J ftT (A- 1 |x a -x^|), (83) 



Ae 2 X 2 S , 

with V the volume of the system, and S the area perpendicular to the applied field. In 
the limit of low density one only has to account for the self-interaction (a = (3) and the 
nearest-neighbor interaction, so that the free energy can be approximated by 



F = BH C1 



1+ 4^ o(A ~ la) 



(84) 



where we used the standard result that H Cl = ln(fi;)/(4eA 2 ) P7[| . This equation shows 
that by raising the applied field an increasing part of the field energy is used to overcome 
the repulsive interaction between the magnetic vortices. This is represented by the last 
term in fl34|), which not only contains a factor of B, like the first term, but in addition 
depends on the lattice spacing a. For increasing fields a decreases, implying that indeed the 
nearest-neighbor interaction term becomes more important. This observation can be nicely 
illustrated by solving the thermodynamic relation H = dF/dB in terms of the magnetic 
induction B, which is related to the vortex density via B = {^/e)p^ (see Fig. |j). For fields 
slightly larger than H C1 , there is a rapid increase in the vortex density which for larger fields, 
when the interaction between vortices becomes important, is leveled off. 

This indicates that the existence of a mixed state is due to the repulsive magnetic inter- 
action between magnetic vortices. If one takes away the interaction term in the free energy, 
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one finds a magnetization curve resembling that of type-I superconductors, which do not 
support a mixed state. The picture is confirmed by the fact that the two critical field val- 
ues H Cl and H C2 , marking the boundaries of the mixed state, become closer to each other 
the smaller the value of the Ginzburg-Landau parameter k — A/£ is. With £ kept fixed, a 
smaller k value implies a smaller penetration depth which results in a stronger shielding of 
the repulsive magnetic force. 

There exist dirty superconductors which have a Ginzburg-Landau parameter close to the 
critical value k = l/y/2 separating a type-II from a type-I superconductor. Such a value 
indicates that the screening of the two forces which are present outside the London limit, are 
of the same order of magnitude. These so-called type-II/1 materials exhibit a remarkable 
experimental phenomenon. If the lower critical field H Cl is crossed from below, a whole flux 
lattice jumps in instead of single vortices being produced one by one as is the case of a deep 
type-II superconductor. This is related to the fact that for these materials the dominant 
force at large distances is the attractive one mediated by the superconducting order field 
28| , while at short distances the repulsive magnetic interaction dominates (see Fig. |^). The 



lattice spacing in the flux lattice produced in such a type-II/1 superconductor is namely an 
equilibrium spacing resulting from the balance between the repulsive and attractive force. 
It is indicated by a in Fig. ||. 

Type-I superconductors, for which the two critical fields H Cl and H C2 have the same value 
H c , and for which k < 1/ V2, do not support a mixed state because the screening length of the 
repulsive magnetic force for those materials is smaller than the coherence length £ defining 
the core radius. That is, effectively the magnetic repulsion is screened to zero outside the 
vortex core, so that vortices in a type-I material experience only an attractive force. They 
will consequently coalesce, with the result that the entire sample becomes normal. 

The other limit we briefly discuss is that of high- vortex density, where the lattice spacing 
a is much smaller than the penetration depth. Because of the periodicity of the lattice, one 
can consider a single unit cell and the Fourier transform becomes a series 
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/i(x) = 5>(k,) e lk *' x , (85) 
{i} 

where 1 stands for the index pair 1 = (Z 1; Z 2 ), and k/ is the reciprocal lattice vector 
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~, - / _-(Z 1 E 1 + / 2 E 2 ) ) (86) 

with E 1 = (-^, — |) and E 2 = (0, 1) two vectors spanning the reciprocal lattice. The Fourier 
components of the plastic field B ext (x) are 

IT xt (kO = -^e ik ^, (87) 

with x a the position of the vortex we are considering. This leads to the free energy 



F — - 
2 



1 /7T\ 2 „ _ 1 



In deriving this we used the fact that all vortices give the same contribution, so that it 
suffice to consider a single vortex, say, located at the origin x a = 0. The term with 1 = in 
( |55D represents the homogeneous part. The remaining sum is diverging for large 1, but this 
is artificial, stemming from the fact that we are in the London limit where the core radius 
is taken to be infinitesimal small. As before, we take an ultraviolet cutoff of the order of the 
inverse coherence length \m^\ = l/£. Since in the limit of large vortex density A/a >> 1, 
one can for 1^0 neglect 1 compared to A 2 k 2 in the denominator of (|88|). In this way, 
approximating the sum by an integral, one finds for the free energy 

„ B 2 1 /tt\ 2 r d 2 k 1 



2 2 Ve/ 7fe min <|k|<fc max (27r) 2 (Ak) 2 

= T + 2^ ln (H' (89) 

The infrared cutoff is taken of the order of the inverse lattice spacing fc m ; n ~ I /a. To be 
precise, 



■k-r' , r p ^H xim6 H' (90) 



where (3 ~ 0.381 is a numerical factor which is adjusted so that the approximate 



calculation is in accordance with the exact result. One thus finds for the free energy 
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B 2 1 f B \ 

F m -BH C2 In , (91) 

2 An 2 2 \H C2 J' 

which in the limit B — > H C2 reduces to the free energy of the normal phase, as it should. 

Intuitively, the phase transition at H C2 may be pictured as follows. When increasing 
the external field, the Abrikosov flux lattice becomes denser, and the spherical cross-section 
of the normal-conducting cores will be continuously deformed into a hexagonal one. This 
picture is nicely confirmed by numerical calculations (see Fig. ||). Precisely at H C2 the 
magnetic vortices are densely packed, and the superconducting cell borders are squeezed to 
zero thickness. In this way the mixed state becomes a homogeneous, normal conducting 
state. 



VII. THE NORMAL PHASE 

In Sec. we saw that the normal conducting phase is described by a disorder theory, 
consisting of a |?/>| 4 theory in the broken-symmetry phase. This theory possesses again 
topological defects, viz., vortices, known from the theory of superfluid 4 He In this 

section a dual formulation of the \ip\ A theory is considered in which the grand canonical 
ensemble of closed vortices is described by a field theory. The theory turns out to be the 
original Ginzburg-Landau model (in the high-temperature phase). This is not surprising 
because a dual map of a dual theory should recover the original model. 

We proceed in the same manner as in Sec. |V| and consider a vortex ending at the point z 
inside the system. To describe such an object we would introduce in the London limit (|78|) 
of the \ip\ 4 theory the operator 

W(L Z ) =e- 1 2 v2 f d H° P fe Xp (l v i J cPxVj.O*^, (92) 



where the first factor is a singular contribution due to the core of the vortex, cf. (pRf). In 
(|92|), the factor n/e stems from our normalization choice of the phase of the ip field, and 
the vector field P , first introduced in the context of superfluid 4 He M, should be such that 
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V x O p yields a delta function along the line L z which starts at the point z and runs to 
infinity along the trajectory y(s): 



with the divergence of this equation yielding a delta function at the point z. But this is 
impossible since for a regular field V- (V x 6 P ) = 0. Hence, the only regular 6 P field that can 
be constructed is one representing a closed, or infinitely long vortex L. To understand this 
in another way ||, imagine a sphere E surrounding the hypothetical endpoint of a vortex 
at z, with a little hole <9E where the vortex pierces the surface. While on the one hand the 
loop integral § dTi dxidtf gives the vortex strength 2tt, the surface integral / s cPxieijkdidjj, 
on the other hand, vanishes since 7 is regular on E. This contradicts Stokes' theorem which 
states that both integrals should be equal and leads to the conclusion that vortices in a |t/>| 4 
theory can only form finite closed loops, or infinite loops which are so to speak "closed at 
infinity" . 

It should be noted that the above argument does not apply to the Ginzburg-Landau 
theory. A magnetic vortex could start in a given point by simply introducing a magnetic 
monopole there. When described with help of a gauge potential, the monopole is inevitably 
accompanied by a Dirac string. Choosing the Dirac string to pierce the surface E surrounding 
the monopole at <9E, one obtains § dTj dxiAi = / E d 2 Xi€ijkdjAi s , where the right-hand side 
measures the magnetic flux through the surface E (see Fig. |7|). Both sides of the equation 
yield the same result, in accordance with Stokes' theorem. 

For a closed loop L, the expectation value of W(L) is given by 



J [ 2 J Ve ' J J 

The integration over 7 can be carried out by substituting the field equation of the Goldstone 
field 




(93) 




(94) 




(95) 
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where A is the correlation function (|53|) with vha = 0. This yields an expression 

(W(L)) = exp |-^ 2 J cPx J d 3 y [V x p (x)], A (x - y) [V x p (y)],} , (96) 

very similar to the one obtained for a magnetic vortex in the superconducting phase, Eq. 
(0). Using (H), we obtain, cf. (0) 

(W(L)) = e" Mlv|i| , (97) 

with \L\ the vortex length and, cf. (|6ll) 



2 'Kr 



7Tf 

"2 ""V V 



M w = — ln(J-f- (98) 



its energy per unit length. In analogy with the previous calculation, we have taken the mass 
\m^\ of the scalar field appearing in the theory as ultraviolet cutoff. Since 7 is massless, 
My/ diverges in the infrared. This is regularized by introducing a small mass \x. 

In order to find the dual theory which gives a field theoretic description of the vortex loop 
gas, we rewrite the expectation value (W(L)) via a Hubbard- Stratonovich transformation 
as 

1. 



(W(L)} = J £> 7 £>bexp{-y d 3 x 



V + ivb ■ - V7 - e l 

2 \e 



(99) 



The integral over the 7 fluctuations now yields the constraint V ■ b = 0, demanding b to be 
the rotation of a vector field, b = V x A. This gives 

(W{L)) = JvAexpi^- J d 3 x -(V x A) 2 -2mvA- J J, (100) 

where Jj(x) := § L dy i 5{'s. — y) describes the closed vortex. It is natural to interpret the 
fluctuating massless gauge field A as the electromagnetic gauge field. This identification 
yields the relation we alluded to above 

2nv = 2e (101) 

between the expectation value v of the disorder field ip and the electric charge 2e. It is the 
analog of relation ( |T9"D between the expectation value w of the superconducting order field 
(f) and the coupling constant g = {Ti/e)mA of the dual theory: 
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2Tiw = g. (102) 

The expectation value (|100| ) we now recognize as the Wilson loop. Since v vanishes as 
T approaches the critical temperature from above, the coupling constant vanishes at the 
critical point and A decouples from the vortex. Precisely the same phenomenon happened 
with the magnetic vortex in the superconducting phase. Adding a gauge-fixing term to fllOOp 
and carrying out the integration over the gauge field, we obtain 



(W(L)) = (e 2 ^^^) =exp [-2e 2 | d 3 xd 3 y J,(x)Aj(x - y) J,(y) 



(103) 



where A°- is the correlation function ( |51~D with vtla = 0. Because we consider closed vortices, 
for which V • J = 0, only the first term of the correlation function contributes, so that the 
result is independent of gauge choice and is given by the previous expression fl96Q. 

If we consider a loop gas of vortices, we recover, following the same line of arguments as 
in Sec. [V], the (normal phase of the) original Ginzburg-Landau model (|1|). As expected, the 
dual map of the dual theory gives back the original model. 

Since the expectation value v of the ip field vanishes when the critical point is approached 
from above, it follows that also the energy ( PSQ tends to zero here. This supports the picture 
that the phase transition in a pure |?/>| 4 theory is associated with the proliferation of vortices 

i- 

The question arises now How can we physically understand a vortex loop in the normal 
phase? As is well known from 4 He physics [pTU l, inside the core of a 4 He vortex, the super- 
fluid order parameter vanishes, indicating that the core consists of normal fluid. Translated 
into the present context, where a finite expectation value of the disorder ip field indicates 
the onset of the normal conducting phase, we are to interpret its vanishing inside the vortex 
core as the absence of the normal phase, i.e., as the presence of the superconducting phase. 
If one takes the Maxwell equation 

V x h = 2ej (104) 
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as the defining equation of the electromagnetic current j, one has to view the purely imag- 
inary object 2niv J in ( |100| ) as describing an electric circuit. (The analog between vortices 
and electric currents was first pointed out by von Helmholtz.) At first glance the presence 
of the factor i seems strange. In particular, the Biot-Savart law for these currents yields the 
opposite sign from what is usually the case: we find two parallel currents repelling instead 
of attracting each other because they carry an imaginary charge. (This can be checked by 
considering the interaction energy E- m t between two parallel vortices 

E int = 2ttV J dxi J d y< -_ — - (105) 

which is positive.) But a closer inspection reveals that this has to be the case. Remember 
that an electric circuit generates a magnetic moment 

m = eJ d 3 x(x x j) (106) 

orthogonal to the surface enclosed by the loop. (Note that we defined the current without 
a charge factor — 2e in our case — included, that is why the equation for m contains a factor 
|2e = e). If we take two real-life circuits I and II as sketched in Fig. |8|, where the lower 
laying loop I is held fixed while the upper one is free to rotate, then the latter would settle 
such that the two magnetic moments point in the same direction. A state where these 
real-life electric circuits are condensed would inevitably be connected with a permanent 
magnetization. Due to the fact that vortex loops carry an imaginary charge, there is an 
antiferromagnetic coupling rather than a ferromagnetic one between the loops, and a vortex 
condensate has zero magnetization. 

We note that because of the imaginary charge, the local field generated by a vortex loop 
L is also purely imaginary as follows from Ampere's law, 

h(x) =*lfdyx j^j^ = i|vfi(x), (107) 

where Q is the solid angle that the loop subtends at x (see Fig. |9|). The same result can be 
obtained from the dual theory, bearing in mind that the local field, apart from a factor i, 
can be identified with the gradient of the phase variable 7 
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h = -iVj, (108) 

see flTT|). Rewriting the field equation fl9"5D for 7, we find that this field can be related to the 
solid angle in the following way M 

7 (x) = f / d\ { ^^l = _ J_ n> (109) 
2tt J s |x — y| d 2tt 

where d 2 yi is an element of the surface S spanned by the loop. Together with ( |108| ) this yields 
the previous result Ql07| ). [The magnetic moment density, or magnetization, is represented 
in the dual theory by 6 P . This follows from the fact that according to ( PH) a closed vortex 
couples to the magnetic field V7 via 6 P .] The order parameter V r (x) of the normal state 
essentially measures the angle Q 

K(x) = e l7r7(x)/e = e iQ(x)/2 . (110) 

We recall that the operator V(x) was constructed by putting a magnetic monopole at x. 
With this kept in mind, Eq. ( |110| ) becomes obvious: Q/2 is the magnetic flux through the 
closed vortex L emanated by the monopole. As a last remark we note that one can chose 
two topologically different surfaces spanning the loop L (see Fig. [Tt]). Both lead, however, 
to the same phase factor because Q differs only by a factor of Air. 

VIII. SUPERCONDUCTING ORDER PARAMETER 

In the previous section it was argued that the transition to the superconducting phase 
could be understood as a proliferation of (closed) vortices of the pure theory. This 
was concluded from the behavior of the real-space representation fl94|) of the single loop 
operator W(L), which was shown to develop an expectation value when T approaches the 
critical temperature from above. The question naturally arises Is W(L) related to the 
superconducting order field 0? To answer this question we have to investigate how W(L) 
is described in terms of the variables of the Ginzburg-Landau theory. To this end we study 
the object 
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0{L Z ) = e ie(z) e 2ie S ^a(x)-e p (x) ; 



111) 



where 6 is the phase of 0, the plastic field E p describes a static charge 2e at z, 



V-E p (x) = <5(x-z) 



and L z is a line emanating from z and running to infinity, 



(112) 



Ef= I d yi 6(x-y). 



(113) 



The second factor in ( 111 ) is incorporated in order to render the operator gauge invariant. 
Indeed, under a gauge transformation 



A(x) -> A(x) + VA(x), 0(z) -> 6{z) + 2eA(z), 



114) 



so that 



0(L.) ^0(L z )exp 



2ieA(z) + 2ze ^ d 3 x VA(x) • E p (x) 



0(L Z ), 



(115) 



where in the last step we performed an integration by parts. To bring out the gauge invari- 
ance of 0(L Z ) more clearly we write it in the equivalent form 



0(L Z ) = exp 



-i J d 3 x(V6 - 2eA) • E 1 



(116) 



We will be working in the low-temperature phase of the Ginzburg-Landau model, where the 
gauge field is massive. We are interested in the expectation value 



(0(L Z )) = J VAV60(L z )e- H , 



117) 



with H the Ginzburg-Landau Hamiltonian (0) in the London limit. Since both integrations 
are Gaussian, they are easily carried out to yield 



(0(L m ))=exp(~J d 3 xd 3 y 
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n(x) A(x - y) n(y) + (2e) 2 E? (x) A(x - y) E?(y) 



118) 



with n(x) = 5(x — z) the charge density. This expression closely resembles the one we found 
for the operator V in (|57|), describing a magnetic vortex, which can be rewritten as 
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(V(L Z )} = exp j-i J d 3 xd 3 y [i n (x) A(x - y) n(y) + ^Sf (x) A(x - y) 5f(y)] }, (119) 

where Bf is defined so that it contains no factor n/e, Bf = (ir/ejBf . We see that ( |118| ) 
can be obtained from ( |119| ) by simply replacing the high-temperature expectation value v 
of the ip field by the low-temperature expectation value w of the <ft field, and by replacing 
the "magnetic" coupling g by the electric coupling 2e. In this sense the operators O and V 
are dual to each other. 

We continue to discuss the behavior of the expectation value ( |118| ) in the two phases. 
In the high-temperature phase we argued that there can only be closed vortices. This we 
achieve by setting n(x) to zero in ( 118|) , so that only the last term survives. In fact, using 



the relation irv = e, we recover the right-hand side of (pop. That is, the operators O(L) and 
W(L) are the same in the high-temperature phase: 

(0(L)) = {W(L)} = e ~ MwW , (120) 

with \L\ the length of the vortex loop and My/ given in (pip. 

In the low-temperature phase the (electric) vortices are condensed and for that reason 
not existing as physical excitations, only the endpoints are physical. The plastic field E p in 
( |118|) can then be written as a gradient of a potential U p , 

E p = -W p , (121) 

with V 2 £7 p (x) = — n(x). Taking a positive charge 2e at z and a negative one — 2e at z, we 
obtain for the correlation function 



(0(L z )0*(L,))=ex P 



J d 3 xd 3 yn(x)A°(x - y)n(y) 



(122) 



2w 2 

where now n(x) = <5(x — z) — 5(x — z), and A°(x — y) is the massless scalar correlation 
function. For x = y we have again a diverging self-interaction which is irrelevant and can 
be eliminated by defining a renormalized operator O r in the same way as we did before in 
(H). We then find 
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(0l ( z)0 ; W ) = exp(^^). (123) 

This low-temperature expression is completely analogous to the one in the high-temperature 
phase for the correlation function (V T (z)V* (z)) , Eq. (|6^). Using the relation e = irv, we can 
write the latter as 

(V;.(z)V;- W ) = exp(J- 3 ^). (124) 

For large separation (O r (z)0*(z)) — > 1, implying that O r (z) develops an expectation value 
in the superconducting phase. That is, O r (z) is the superconducting order parameter. Being 
gauge invariant this operator makes no statement about the local U(l) symmetry. Referring 
back to the first representation ( |111| ) of the superconducting order parameter, we find that 
a non-zero expectation value indicates that the global U(l) symmetry parameterized by a 
constant transformation parameter A is spontaneously broken. 

In closing this section, we remark that the superconducting order parameter can also be 
represented in the dual theory. The result is that E p appears in (173) in the combination 
V x h — 2eE p with the fluctuating h field. To derive this result it is prudent not to proceed in 
the manner we exploited before to obtain the dual theory (fT4j) and linearize the gauge-field 
fluctuations, but instead linearize the 9 fluctuations. The result can also be inferred using a 
duality argument, remembering that the disorder parameter was incorporated in the original 



theory by the combination V x A — B , see Eq. ([46]). 



IX. DISORDER FIELD THEORY FOR THE SUPERCONDUCTING PHASE 

TRANSITION 

We have argued that the dual description of the Ginzburg-Landau model is one in terms 
of physical variables, and that it possesses no local gauge symmetry. A disorder field ip was 
identified which vanishes in the superconducting phase, and which develops an expectation 
value in the normal conducting phase, thereby breaking a global U(l) symmetry. We recall 
that central to Landau's theory of continuous phase transitions is the presence of an order 
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parameter which signals through its vacuum expectation value whether or not a certain 
symmetry is broken. This is precisely what the disorder field ip does. On that ground ip 
is ideally suited to formulate a Landau type of description of the superconducting phase 
transition P,|18|j. To understand why the Ginzburg-Landau theory itself is not well suited to 



do this, we note that it has a local gauge symmetry. According to Elitzur's theorem [ IIJ such 
a symmetry can never be broken, so that for a local symmetry no order parameter exists 
in the sense of Landau. This may be one of the reasons for the fact that no infrared stable 
fixed point was found within the Ginzburg-Landau formulation of the superconducting phase 
transition although it is generally accepted that the transition is of second order in 



the type-II regime, and thus should possess such a point. 

Below the transition temperature we saw that the dual theory consists of a |^| 4 theory 
coupled to a massive vector field. Above T c this field decouples, and the disorder field %fj 
simultaneously develops a vacuum expectation value. Despite the apparent difference in the 
description of the low- and high-temperature phase, it is readily argued that at the mean- 
field level the critical behavior is governed by a simple |?/>| 4 theory. This can be seen by 
integrating out the massive vector field in the low-temperature phase. Apart from irrelevant 
terms, this leads to only a change in the coefficients of the |^| 4 theory, no additional relevant 
terms such as are generated. Explicitly, 



H^s = I d x 



1 2 , / 2 2 171 A \ I / 1 2 i | 9 4 



(125) 



In deriving this effective Hamiltonian we used dimensional regularization; (irrelevant) higher- 
order terms were omitted. We note that all contributions stemming from the vector field 
h vanish in the limit T approaching T c from below, so that f^, e ff reduces to ([77]) in this 
limit. It is well known that a \tp\ 4 theory with a positive coupling has a non-trivial infrared 
stable fixed point and undergoes a second-order phase transition. We therefore conclude 
that at the mean-field level also the superconducting phase transition is of second order if 
the system is sufficiently deep in the type-II regime. The same conclusion can be reached 
starting from the Ginzburg-Landau formulation. 
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Below we will apply renormalization group theory to see if this conclusion holds also 
beyond mean-field theory. Halperin, Lubensky, and Ma ]32] performed this study within the 
Ginzburg-Landau theory. Using an e expansion, they showed that at the one-loop level the 
theory no longer possesses an infrared stable fixed point. They interpreted this as indicating 
that the transition is of first order. This conclusion was in accordance with results obtained 



by Coleman and Weinberg p4| who studied the electrodynamics of massless scalar mesons 
in four dimensions and discovered that at the one-loop level the photon acquires a mass. A 
study of the effective action shows a precocious onset of the Higgs mechanism with a sudden 
appearance of a finite photon mass. This is typical for a first-order transition. Only by 
artificially enlarging the number of components of the complex scalar field appearing in 
the Ginzburg-Landau model did Halperin, Lubensky, and Ma find an infrared stable fixed 
point, provided this number is taken to be larger as 183. However, the corresponding critical 
exponent rj, which determines the anomalous dimension of 0, depends on the gauge-fixing 
choice and is therefore unphysical. This should not come as a surprise since <ft itself is not 
gauge invariant, and therefore not physical. 

We shall carry out the renormalization group theory within the dual formulation which, 
being casted in terms of physical fields, does not suffer from the flaws of local gauge invari- 
ance fl[]. The dual theory (|[4|) involving a massive vector field is perturbatively renormal- 
izable in four dimensions (D = 4) ||35|| , so that usual perturbation theory can be applied to 
calculate the critical exponents. However, the derivation of the dual theory (|74|) from the 
Ginzburg-Landau model hinged on the fact that the number of space dimensions is three, 
for which the dual object e^kdjA^ is a vector. In other words, the dual theory describes the 
superconducting phase transition only in three dimensions. For this reason, we carry out the 
renormalization group theory in fixed (D = 3) dimension, and not in D = 4 — e dimensions 
as is often done. The fixed-dimension approach to critical phenomena was introduced by 
Parisi, who applied it to a pure \ip\ A theory p6|JT4]] . The method makes explicitly use of 
the fact that near the critical point the system has only one relevant length scale, viz., the 
correlation length which diverges at this point. This length is used to convert dimensionful 
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coupling constants into dimensionless ones. 

In the present setting the relevant scale is the (renormalized) inverse mass m^ 1 . (The bare 
mass vanishes as m\ ~ T c — T at the critical temperature T c ). We know from the Ginzburg- 
Landau theory that the bare penetration depth also diverges at T c , viz., Ao ~ (T c — T) -1 / 2 . 
However, the renormalized length A should not constitute an independent diverging length 
scale. We will see below that this is indeed the case. As usual, the critical exponents are 
computed in the symmetric phase of the model, which in the present context corresponds 
to the superconducting phase. 



We write the bare Hamiltonian ( 74]) as a sum of the renormalized Hamiltonian and 
counterterms 5H 



5H = / d 3 x 



[Z^ - 1)|(V - igh)^\ 2 + (V4 j0 - mJ)|Vf + u(Z u - 1)|W 4 



+i(Z h -l)(Vxh) 2 + l(Z h m 2 4)0 -m 2 4 )h 2 



(126) 

(All quantities appearing in (f73|) should have been given an index to indicate that they 
refer to bare quantities. For convenience we rescaled the h field by a factor m^: h — ► m^h.) 
The renormalized objects are related to the bare ones via 

hi = Z- l/ \„ g = Z g l Z,X /2 9o, ^ = Z~ l/2 ^ u = Z- 1 Z 2 4) u . (127) 

It is straightforward to calculate the one-loop diagrams. The correlation functions can be 
read off from the Hamiltonian (0). With a wiggly line denoting the correlation function of 
the vector field, and a straight line denoting the one of the ip field, it follows that 

1 ( h - W) ( i28a ) 



k 2 + m\ 



k 2 + m\ ' 



;i28b) 



The correlation function ( 128a|) reflects the constraint V • h = which we imposed upon 



the fluctuating field h. The correlation function has been obtained by representing the S 
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function <5(V • h) as exp[— (V • (h) 2 /(2£)] with ( taken to zero at the end. We find for the 
diagrams depicted below 

2 9 2 l2 



3n m^p + rriA 



(129a) 




2 q 3 

i- h (129b) 

Sir m^p + rriA 



1 " 2 -k 2 fc-^V (129c) 



24^/7^ V J k 2 



We considered these particular diagrams to illustrate the following observations. The first 
diagram yields a Z^p factor given by 

2 a 2 

Z i> = l + y - . (130) 

From the second diagram we extract the factor Z g , 

z 9 = i + ^- — x — ' ( 131 ) 

67V m^p + VTiA 

which turns out to be equal to Z$. To this order, the renormalized coupling constant g is 
thus related to the bare one simply via 

9 = zl /2 go. (132) 

This means that the minimal coupling to the vector field is preserved at the one-loop level. 
In fact, it is preserved to any order in the loop expansion. The reason is that the Ward 
identity which guarantees that the minimal coupling is always maintained in the case of a 
massless vector field, also operates in the massive case ||35|| . 

The diagrams ( p.29cj ) show that the one-loop contributions to the self-energy of the vector 
field is transverse. This too remains true to all orders in perturbation theory thanks to the 
Ward identity. The mass term of the vector field is consequently not renormalized and does 



1/2 

not need a counterterm. That is, rriA = Z h mA,o, so that the last term in ( |126|) vanishes. 
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It follows that the critical exponent 7h is unaffected by the fluctuations, and retains its 
mean-field value 7^ = 1. 

Incidently, the electric charge does not renormalize in the dual theory since g = (7r/e)m^, 
and both g and rriA renormalize in the same manner. 

We now come to an important observation related to the fact that in the dual theory 
the inverse penetration depth plays the role of a mass as well as of a coupling constant 
because g = (7i/e)mA, with e a constant. The standard definition of the critical exponent v 
which determines how the correlation length m^ 1 diverges when the temperature approaches 
T c : ~ (T c - T)~ u , is 

1 dln{mi ) 



v <91n(m</,) 
In our case this can be rewritten as follows 



(133) 



= ^ (134) 
because toL q ~ m AQ near T c . The (3 function is defined by the equation 

^f);= m ®l- , (135) 

with the properly scaled coupling constant g 2 := g 2 /m^. By virtue of ( |134j ), this can be cast 
in the form 

m=f (. 1 + l + lki ^ ) a J^Ay (13 6) 

where ■jh(g 2 ,u), with u := u/m^, is the function 

d 

-fh(g 2 ,u) := m A — -ln(Z h )| MOigo (137) 

which yields the critical exponent rjh when evaluated at the critical point. Without the 
explicit mass dependence, the coefficient of the g 2 term in the f3(g 2 ) function would be —1, 
implying that the origin g = is an ultraviolet stable fixed point. In ( |136[ ), however, the 



coefficient is —1 + l/v which is positive if v < 1. In this case, the origin g 2 = becomes 
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infrared stable and the coupled theory reduces to a pure |?/>| 4 theory. The best estimate for 
v available from summed perturbation theory at fixed D = 3 |37| gives v « .6695, which is 
smaller than one. Hence, the trivial fixed point g* = is infrared stable (see Fig. [11]). This 
situation differs dramatically from that in the Ginzburg-Landau theory where the coupling 
e to the electromagnetic gauge field has an infrared stable fixed point away from the origin, 
see Fig. [12], and the corresponding value e* is too large to allow the coupled system to 
develop an infrared stable fixed point. 

To recapitulate, the dual theory of the superconducting phase transition possesses an 

2 

infrared stable fixed point given by g* =0 and u* = u WF , where iiyjF * s the Wilson-Fisher 
fixed point of a pure 4 theory with reversed temperature axis. The critical exponents of 
the ip field are the ones of a superfluid. The critical exponents pertaining to the h field, 
which physically represents the local induction, have their mean-field values. In particular, 
Uh = 1/2. This exponent reveals that the magnetic penetration depth diverges near T c 
as (T c — T) -1 / 2 , meaning that inside the critical region the empirical formula A ~ [1 — 
(T/Tc) 4 ]- 1 / 2 found outside this region remains unchanged. 

2 

A last point of interest is the Gaussian fixed point, corresponding to g* = 0, ii* = 0. This 
fixed point is infrared stable in the g 2 direction, and unstable in the u direction. It describes 
a loop gas of free vortices. At this point the phase transition changes from second to first 
order, i.e., it is a tricritical point, the existence and location of which was first established 



in Ref. ||38|| . At the level we are working, the critical exponents characterizing this point are 
Gaussian. A \ip\ 6 term which should now be included will generate logarithmic corrections. 
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APPENDIX: 



Here we give a lattice version of the arguments presented in Sec. [HI] , showing the equiva- 
lence of a complex field theory and a loop gas [|,nj. We start with the free Hamiltonian (|2T| 



defined on a hypercubic lattice in D space dimensions. (We use the same lattice notation 
as in Sec. |TV| ). The corresponding energy expression we write as 

£ = ]Ta D -V(x)A(x,y)^(y), (Al) 

where 

A(x, y) = (2D + mV)5 x , y - ]T (<5 x , y+ai + <5 x , y _ ai ) , (A2) 

i 

or in matrix notation 

A = (2D + m 2 a 2 )I - 5. (A3) 

Here, I is the identity operator and S is the so-called step operator. The matrix element 
S*(x, y) is 1 if the two lattice sites x and y are neighbors, and zero otherwise. From (ST) 
we immediately read off the lattice correlation function 

G(x) =a 2 - D A- 1 (0,x). (A4) 

The operator A -1 we expand in a von Neumann series 

where we introduced the abbreviation 



1 m 2 a 2 , . . 

° = 1 + -2D- < A6 » 

From the definition of the step operator S it follows that S N (x,y) is equal to the number 
of paths which go from site x to site y in N steps. By virtue of this we can write 

oo 

A- 1 (0,x) = ^-E^( x )^ (A7) 
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with Pjv(x) the probability defined by the number of paths that go from to x in N steps 
divided by the ones that just start in and contain N steps. This last number is equal to 
(2D) N . The probability P/v(x) plays a central role in the theory of random walks. A little 
thought reveals that it satisfies the recurrence relation 



Pv+iM 



— ^Pv(x + ai), 



(A8) 



with the initial condition 



Pnfx) = 5, 



0,x- 



(A9) 



To solve ( |A8|) we introduce the Fourier transform: 



D I a ft ikx r> /i \ 

a / — x Pv(k). 



r/a (2n) D 

The recurrence relation then yields for P/v(k) 



Pv+i(k) = — ^cos(/cia)Pjv(k), 



(A10) 



(All) 



with Po(k) = 1 as follows from the initial condition (|A9|). In this way we obtain as solution 



Pv(x) = a 



D 



*/* d D k 



ikx 



N 



(A12) 



We are now in a position to take the continuum limit (a — > 0). Let us first concentrate 
on the probability (A12) and assume that the length s of a path on a lattice is measured in 
some — as yet arbitrary — unit b. That is, a path of iV steps has length s = Nb. Since in the 
limit a — ► 



D 



^COs(/CjC 



N 



s/b 



cxp 



2Db 



(A13) 



we have to take the ratio a 2 /b fixed in order to obtain a non-trivial continuum limit; we will 
take a 2 /b = 2D. In other words, b ~ a 2 in the scaling limit. In this limit the probability 
density p(x, s) := lim a _^ P N (x.)/a D becomes 

d D k 



P( x > s ) = J 



{2tt) d 
46 



e ik-x g -sk 2 



(A14) 



so that the correlation function ( \K^) , with A 1 given by ( |A5|) , can be written as 

a 2 a 00 1 / m 2 a 2 \~ s,b 

G(x) = lima-«A-'(0, x) = Jim i£ £ -Lfl,(x) 1 ' " ' 



a-o v ' ' a^o 2D ^ a D y ' \ 2D 



co 2 







e- sm p(^s) } (A15) 
where we used the fact that a — > 1 and 



a 2 00 



00 



2D N=Q Jo 

when a —>■ 0. But with p(x, s) given in (|A14| ), this is precisely the Schwinger representation 
(f23|) of the correlation function. 

A similar derivation can be given for the partition function. To this end we consider the 
probability P^ defined by the number of closed oriented paths containing A" steps divided 
by (2D) N . It is related to the probability -Pat(x) previously introduced via 

P N (0) = NP C N , (A17) 

where the factor N arises from the fact that one can start traversing a given oriented loop 
at any of the A" lattice sites visited by the loop. It then easily follows that 

V tJ<a N - / -e~ sm2 / f-4re- sk2 , (A18) 
^ a D Jo s J {2tt) d v ; 

which is the right-hand side of the first equation in (p7|). This leads to the same repre- 
sentation of the partition function we found before with help of Schwinger's proper-time 
representation. The steric repulsion can also be included on the lattice, but this will not be 
done here. 
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FIGURES 

FIG. 1. Magnetic field (thin lines) of a positively charged monopole. The magnetic flux is 
provided by the Dirac string (thick line) located along the negative z axis. 

FIG. 2. Biot-Savart interaction (wiggly line) between two line elements dxi and dyi of a mag- 
netic vortex (straight line). 

FIG. 3. Triangular Abrikosov flux lattice. The two vectors vi and V2 span the unit cell, 
containing one flux tube. 

FIG. 4. Magnetic induction B versus applied magnetic field H in arbitrary units for a type-II 
superconductor. 

FIG. 5. Schematic representation of the interaction potential V(r) between two parallel mag- 
netic vortices for a type-II/1 superconductor. 

FIG. 6. Contour plots of the local induction h in an Abrikosov flux lattice for a type-II super- 
conductor with k = 1.3. The first plot is at B/H C2 = 0.1 and the second is at B/H C2 = 0.9 (after 
Ref. ||). 

FIG. 7. A magnetic monopole with its string piercing the surface £ at <9£. 

FIG. 8. Two magnetic moments mj and 1117-7- interacting with each other. The first one is kept 
fixed, while the second is free to rotate. 

FIG. 9. Solid angle f2 that the vortex loop J subtends at the point x. 

FIG. 10. Two different surfaces spanning the loop L. The flux through the surfaces differ by a 
factor 47T. 

FIG. 11. The j3 function for the dual coupling constant (squared) g 2 = g 2 /m^, showing that 
the origin is an infrared stable fixed point. 
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FIG. 12. The f3 function for the electric charge (squared) e = g/m^ as obtained from the 
Ginzburg-Landau theory, showing that the origin is an infrared unstable fixed point. 
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FIG. 1 

r-r- dXj 

FIG. 2 
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FIG. 5 




FIG. 6 
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FIG. 8 
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Q(x) 



FIG. 9 
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FIG. 10 
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